Domain wall motion of magnetic nanowires under a static field 
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The propagation of a head-to-head magnetic domain-wall (DW) or a tail-to-tail DW in a magnetic 
nanowire under a static field along the wire axis is studied. Relationship between the DW velocity 
and DW structure is obtained from the energy consideration. The role of the energy dissipation in 
the field-driven DW motion is clarified. Namely, a field can only drive a domain-wall propagating 
along the field direction through the mediation of a damping. Without the damping, DW cannot 
propagate along the wire. Contrary to the common wisdom, DW velocity is, in general, proportional 
to the energy dissipation rate, and one needs to find a way to enhance the energy dissipation in order 
to increase the propagation speed. The theory provides also a nature explanation of the wire-width 
dependence of the DW velocity and velocity oscillation beyond Walker breakdown field. 

PACS numbers: 75.60.Jk, 75.60.Ch, 85.70.Kh, 74.25. Ha 
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Manipulation of domain wall (DW) of magnetic 
nanowires by a field 0,i and/or by a current 

0, H, S E3] has recently attracted much attention because 
of its academic interest and potential applications in in- 
formation storage llj . Many interesting phenomena were 
discovered with limited understandings. For a tail-to-tail 
DW or a head-to-head DW as shown in Fig. 1 in a mag- 
netic nanowire with its easy-axis along the wire axis, the 
DW will propagate in the wire under an external mag- 
netic field parallel to the wire axis or a current through 
the wire. The propagating speed v of the DW depends on 
the field strength [1,13] and/or the current density pH Hoj. 
For the field-driven DW motion, there exists a so-called 
Walker's breakdown field Hur fl2l| . The DW velocity v 
is proportional to the external field H for H < H\,y 
and H ^> Hw- DW velocity decreases as the field in- 
creases between the two linear H-dependencc regimes. 
The linear regimes are characterized by the DW mobility 
fx = v/H. For H ^> H\y, the DW velocity may also os- 
cillate with timepjl. Experiments also found that the 
DW velocity is sensitive to both DW structures and wire 
width0,|3,|3|. Magnetization dynamics is governed by a 
nonlinear differential equation, called Landau-Lifschitz- 
Gilbcrt (LLG) equation. Spatial homogeneous solutions 
of this equation related to the optimal magnetization re- 
versal of Stoner particles may be solved analytically [3], 
but general solutions for the DW motion are not known. 
Our current theoretical knowledge about DW motion is 
either from the seminar work of Walker (T^j or from the 
Slonczewski's one-dimensional modelfTo| in which a DW 
is treated as a rigid body characterized by DW position 
and a cant angle from the demagnetic field of a film. 
Both theories are developed for one-dimensional trans- 
verse DW of a uniaxial magnetic wire with a demagnetic 
field. 

Though much progress have been made on both experi- 
mental and theoretical aspects, DW motion is still poorly 
understood. Slonczewski model is a great simplification 
of DW motion from the LLG equation. This model does 
not contain the detail structure of DW with DW width 
as an input parameter. Thus the theory can explain nci- 
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FIG. 1: Schematic diagram of a head-to-head domain wall 
(DW) in a magnetic nanowire. There are three phases in the 
wire, two magnetic domains I and II and DW of width A. 
Domain I has its magnetization pointing to +z-direction (9 = 
0) and the magnetization of Domain II is along -z-direction 
(9 — 7r). Ill is the DW region whose magnetization structure 
could be very complicate. An external field is assumed to be 
along +z-direction. DW propagates along the field direction. 



ther DW velocity oscillation nor many other interesting 
issues. Walker's theory is based on a particular solution 
of the LLG equation in a special case that does not dis- 
tinguish different types of DWs. It predicts that the DW 
mobility is proportional to the DW width and inversely 
proportional to the damping constant, /i = 7A/0J, where 
7 is gyromagnetic ratio, A is the DW width and a is 
damping constant. Walker's analytical analysis is useful 
in understanding some aspects of experimental findings, 
but it does not provide deep insight into DW motion such 
as the origin of the motion. In fact, Walker's expression 
of the DW mobility is strange because it is singular in 
a. It docs not tell us what will happen when damping 
goes to zero (the limit of lim^^o lim a _>o docs not exist, 
and the limit of lim Q _>o lim#=^# w ^o depends on j3 that 
is a numeric number smaller than 1). Although Walker's 
solution is only correct for a strictly one-dimensional uni- 
axial wire with a demagnetic field and there is no proof 
of its correctness for other types of wires, it was widely 
used to understand DW motion. Thus, it is important 
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to prove whether /i = 7 A /a is true for an arbitrary mag- 
netic wire and to know how the result should be modified 
if it is not. 

In this letter, we shall examine the field-induced DW 
propagation from an angle different from those of Walker 
and Slonczewski. Firstly, we shall, from a simple model 
calculation, show that no static head-to-head DWs or 
tail-to-tail DWs are allowed in the presence of an exter- 
nal magnetic field along the wire axis. Secondly, energy 
conservation principle prevents a DW from propagating 
in the absence of a damping even under a static field, 
and it leads to a general relationship between the DW 
velocity and DW structure. Contrary to the common 
wisdom, DW velocity is proportional to the energy dis- 
sipation rate. Thus, the intrinsic critical field for DW 
propagation in a wire is zero, and the DW velocity is in 
general proportional to damping constant. Finally, our 
theory can reproduce the results of both Walker's and 
Slonczewski's models. Furthermore, the present theory 
attributes a DW velocity oscillation for H 3> Hw to the 
periodic motion of the DW, either the precession of the 
DW or oscillation of the DW width. 

In a magnetic material, magnetic domains are formed 
in order to minimize the stray field energy A DW that 
separates two domains is defined by the balance between 
the exchange energy and the magnetic anisotropy energy. 
The stray field plays little role in a DW structure. To de- 
scribe a head-to-head DW in a magnetic nanowire, let us 
consider a wire with its easy-axis along the wire axis (the 
shape anisotropy dominates other magnetic anisotropics 
and makes the easy-axis along the wire when the wire is 
small enough) which is chosen as the z-axis as illustrated 
in Fig. 1. Since the magnitude of the magnetization M 
does not change in the LLG equation [l4|], the magnetic 
state of the wire can be conveniently described by the 
polar angle 9(x, t) (angle between M and the z-axis) and 
the azimuthal angle <f>(x,t). The magnetization energy 
is mainly from the exchange energy and the magnetic 
anisotropy because the stray field energy is negligible in 
this case. The wire energy can be written in general as 



(1) 



E = J F{9,<p,V6,V<j))d 3 x, 

F = f(8, 4>) + J[( V6>) 2 + cos 2 6»(V0) 2 ] - MH cos 6, 



where / is the energy density due to all kinds of magnetic 
anisotropies which has two equal minima at 9 = and 
n = 0, 0) = f(6 = 7T, </>)), J— term is the exchange 
energy, M is the magnitude of magnetization, and H is 
the external magnetic field along z-axis. In the absence 
of H, a head-to-head static DW that separates 9 = 
domain from 9 = ir domain (Fig. 1) can exist in the wire. 
The domain structure is determined by the equations 



= 2JV 2 6» -% + 2Jsin6»cos6»(V0) 2 = 0, 
69 09 
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In order to show that no intrinsic static head-to-head 
DW is allowed in the presence of an external field (H 7^ 
0), consider a special case that the wire is ID and / is 
a function of 9 only. Then Eq. @ should include the 
external field and becomes 



r d 2 9 df 



- ^ -#Msin(9 + 2Jsin6<cos(9(-^) 2 = 0, 
dz 1 o9 dz 

dd> 

2 J cos 9— = const, 
dz 



(3) 



The second equation requires d<f>/dz = 0. Since <f> = 
const, and 9 = or 7T on the both sides of the DW, 
d<p/dz = 0. Thus, one can multiply the first equation 
by d9 J dz and the integration of the first equation yields 
J(d9/dz) 2 -f{6)- HM cos 9 = C = const. However, for 
a head-to-head DW as shown in Fig. I, C = — /(0) — 
HM in the far left is not compatible with C = —f(ir) + 
HM in the far right. Thus, Eq. §5§ does not have a 
head-to-head DW solution for H ^ 0. In other words, 
a DW in a nanowire under an external field must be 
time dependent that could be either a local motion or a 
propagation along the wire. It should be clear that above 
argument is only true for an intrinsic head-to-head DW, 
and it shall fail with defect pinning that changes Eq. (J3]). 
Indeed static DWs exist in the presence of a weak field 
in reality because of pinning. 

The following results arc based on the fact that a static 
field can be neither an energy source nor an energy sink. 
Without a damping, the system is a Hamiltonian sys- 
tem, and will move on an equal energy contour by energy 
conservation. This viewpoint was important in under- 
standing magnetization reversal of Stoner particles in a 
static fiel dllol , and its extension was used to develop new 
strategies [13 for the magnetization reversal of Stoner 
particles. From the LLG equation [l6|, one can show that 
the energy damping rate is 
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dE_ _ 
~dt ~ 

where H e fj = is the effective field. In regions I and 

II or inside a static DW, M is parallel to H e ff. Thus no 
energy dissipation is possible there. The energy dissipa- 
tion can only occur in the DW region when M is not 
parallel to H e f / . 

The total energy of the wire equals the sum of the en- 
ergies of regions I, II, and III, E = 11 in^i- Ei 
increases while Eji decreases when the DW propagate 
from left to the right along the wire. The net energy 
change Ej + Eji in region I plus II due to the DW prop- 
agation is 



d{Ei + E n ) 
dt 



-2HMvA, 



(5) 
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where v is the DW propagating speed, and A is the cross 
section of the wire. The energy of region III should not 
change much because the DW width A is defined by the 
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balance of exchange energy and magnetic anisotropy, and 
is usually order of 10 ~ lOOnm. A DW cannot absorb 
or release too much energy, and can at most adjust tem- 
porarily energy dissipation rate. In other words, dE d ' t " 
is either zero or fluctuates between positive and nega- 
tive values with zero time-average. Since energy release 
from the magnetic wire should be equal to the energy 
dissipated (to the environment), one has 



2HM vA 



dEru 
dt 



1 + 



M x H P 



(6) 

Thus, the relationship between the DW velocity and the 
DW structure is 



2(1 + a 2 )HA J UI 



rh x H e ff ) d 3 x- 



1 dE 



in 



2HMA dt 



(7) 

where m is the unit vector of M. Eq. (J7]) is our central 
result. Obviously, the right side of this equation is fully 
determined by the DW structure. Time averaged velocity 
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d 3 3 
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where bar denotes time average. 

Because there is no intrinsic static head-to-head DWs 
in the presence of a static field, the first term in the 
right side of Eq. ([7]) will be positive and non-zero since 
a time dependent DW requires m x H e ff ^ 0. Thus 
time-averaged v is always positive. It implies zero intrin- 
sic critical- field for DW propagation which is consistent 
with Walker's result (T3|. A DW can have two possible 
types of motion under an external magnetic field. One is 
that a DW behaves like a rigid body propagating along 
the wire. This is the basic assumption in Slonczewski 
model [3, EH and Walker's solution for H < Hw- Obvi- 
ously, DW energy is time- independent, dE J t " = 0. The 
other one is that the DW structure varies with time while 
it propagates in the wire. It means that the snap-shots 
of the DW at different time are different. In the first 
type of motion, if the DW keep its static structure in the 
absence of the external field, then the first term in the 
right side of Eq. ((7|) shall be proportional to aAAH 2 , 
where a, characterize the DW structure, is a numerical 
number of order of 1 (smaller than 1). This is because the 
effective field due to the exchange energy and magnetic 
anisotropy is parallel to M, and does not contribute to 



the energy dissipation. Thus, in this case, v 



ay A 



H 



with /i 



ay A 
'T+a~? 



To show that a depends on material pa- 



rameters, let us consider the Walker 's[12| model in which 
f = -^M 2 cos 2 6+ ^M 2 sin 2 9 cos 2 0, here K x and K 2 
measure the magnetic anisotropics along the easy and 
hard axes, respectively. From Walker's trial function 



In tan v ' ' = c(t) 
( j>{z,t)=4>{t), 



v(r)dr 



(9) 



where c _1 (i) is the DW width A, Eq. (4} becomes 

f = -^^y [KiM 3 sin 2 </,cos 2 + H 2 M] , (10) 
and DW energy change rate is 
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dt 



dt 



F{0,(j),Ve,V(l))d 3 x^AJA-c(t). (11) 



in 



Substituting Eqs. JTDJ) and dTTJ) into Eq. 0, and to- 
gether with Walker's conjectures of C and <f>, one can eas- 
ily reproduce Walker's DW velocity expression for both 
H < Hw and 3> Hw- For example, for H < Hw = 
aK 2 M/2, c _1 (t) = A = const, and Eq. |J7J) gives 



aA'j 
1 + a 2 



K 2 M sin <f> cos < 
H 



H. 



(12) 



This velocity expression is the same as that of the Slon- 
czewski model[3, [TH in a one-dimensional wire. In 
Walker's analysis, <j> is fixed by K 2 and H through 
K 2 M sine/) cos <j> = —. Using this (f> in the above ve- 
locity expression, Walker's mobility coefficient pi = — 
is recovered. This inverse damping relation is from the 
particular potential landscape in (^-direction. One should 
expect different result if the shape of the potential land- 
scape is changed. Thus, this expression should not be 
used to extract the damping constant [H, 

A DW may be substantially distorted from its static 
structure when H 3> Hw as it was revealed in Walker's 
analysis. Generally speaking, a physical system un- 
der a constant driving force will first try a fixed point 
solution [l8|- It will go to other types of more compli- 
cated solutions only if a fixed point solution is not possi- 
ble. In terms of DW motion, a domain wall will arrange 
itself as much as possible to satisfy Eq. (3). Of course, 
it is known that a head-to-head DW solution is only pos- 
sible for Eq. (2) but not for Eq. (3). Thus, one should 
expect that AI is almost parallel to H e f f except a small 
perpendicular component of order of H. One may as- 
sume M x H e f f ~ bAIH sin ip with b much smaller than 
1. This assumption agrees with the minimum energy 
dissipation principle fl9j since \AI x H e ff\ = MHsitlQ 
when M for H = is used, and any modification of M 
should only make \M x H e ff\ smaller. The DW prop- 
agating speed is then v = b'ajAH/(l + a 2 ), a linear 
relation in both A and H . but a smaller DW mobility 
because of the smallncss of b. This smaller mobility at 
H ^> Hw leads naturally to a negative differential mo- 
bility between H < Hw and H ^> Hw^- In other words, 
this negative differential mobility is due to the transition 
of the DW from a high energy dissipation structure to a 
lower one. This picture tells us that one should try to 
make a DW capable of dissipating as much energy as pos- 
sible if one wants to achieve a high DW velocity. This 
is very different from what people would believe from 
Walker's special mobility formula of inverse proportion 
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of the damping constant. To increase the energy dissipa- 
tion, one may reduce defects and surface roughness other 
than increasing the damping constant. The reason is. by 
minimum energy dissipation principle, that defects are 
extra freedoms to lower \M x H e ff\ because, in the worst 
case, defects will not change \M x H e ff\ when M without 
defects are used. 

In the case that the DW precess around wire axis 
or DW width breathes periodically or both motion are 
present, one should expect both dE J t " and energy dissipa- 
tion rate oscillate with time. According to Eq. (ff|. DW 
velocity will also oscillate. DW velocity should oscillate 
periodically if only one type of DW motion (precession 
or DW breathing) presents, but it could be very irregular 
if both motions are present and the ratio of their peri- 
ods is irrational. Indeed, this oscillation was observed 
in a recent experiment 3 . It should be pointed out that 
the DW precession motion was already detected experi- 
mentally by MOKE broadening technology @. How can 
one understand the wire-width dependence of the DW 
velocity? According to Eq. ([7]), the velocity is a func- 
tional function of DW structure which is very sensitive to 
the wire width. For example, for a very narrow wire, only 
transverse DW is possible while a vortex DW is preferred 
for a very wide wire (large than DW width). Different 
vortices yield different values of \M x H e f/\, which in 
turn results in different DW propagation speed. 

The correctness of our central result Eq. {7} depends 
only on the LLG equation, the general energy expres- 
sion of Eq. (fT|), and the fact that a static magnetic field 
can be neither an energy source nor an energy sink of a 



system. It docs not depend on the details of a DW struc- 
ture as long as the DW propagation is induced by a static 
magnetic field. In this sense, our result is very general 
and robust, and it is applicable to an arbitrary magnetic 
wire. However, it cannot be applied to a time-dependent 
field or the current-induced DW propagation, at least not 
directly. Also, it may be interesting to emphasize that 
there is no inertial in the DW motion within LLG de- 
scription since this equation contains only the first order 
time derivative. Thus, there is no concept of mass in this 
formulation. 

In conclusion, a global view of the field-induced DW 
propagation is provided, and the importance of energy 
dissipation in the DW propagation is revealed. A gen- 
eral relationship between the DW velocity and the DW 
structure is obtained. The result says: no damping, no 
DW propagation along a magnetic wire. It is shown that 
the intrinsic critical field for a head-to-head DW is zero. 
This zero intrinsic critical field is related to the absence 
of a static head-to-head DW in a magnetic field parallel 
to the nanowirc. Thus, a non-zero critical field can only 
come from the pinning of defects or surface roughness. 
The observed negative differential mobility is due to the 
transition of a DW from a high energy dissipation struc- 
ture to a low energy dissipation structure. Furthermore, 
the DW velocity oscillation is attributed to either the 
DW precession around wire axis or from the DW width 
oscillation. 
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